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ABSTRACT 

We have calculated the pressure and energy density in four-flavour QCD using im- 
proved fermion and gauge actions. We observe a strong reduction of finite cut-off 
effects in the high temperature regime, similar to what has been noted before for the 
SU(3) gauge theory. Calculations have been performed on 16^ x 4 and 16^ lattices 
for two values of the quark mass, ma = 0.05 and 0.1. A calculation of the string 
tension at zero temperature yields a critical temperature \fo = 0.407 ±0.010 for 
the smaller quark mass value. 



1 



1 Introduction 

In the pure gauge sector of QCD improved actions have been shown to lead to 
a drastic reduction of systematic errors introduced by the non-zero lattice spacing 
{finite cut-ojf effects) in the calculation of thermodynamic observables ^ 0. This 
is particularly evident in the high temperature regime where analytic calculations 
in the infinite temperature, ideal gas limit show that deviations from the continuum 
result can be drastically reduced already with tree level improved Symanzik actions 
However, also close to and even at an improvement has been observed 
already with tree level improved actions 0. 

At non-zero temperature the finite cut-off effects become visible as a finite size 
effect because the lattice spacing in units of the temperature is given by the temporal 
extent, iV,-, of the lattice, i.e. aT = 1/iV^. In the standard Wilson formulation of 
lattice QCD with staggered fermions these cut-off effects are known to be 0((aT)^ = 
■ In the high temperature, ideal gas limit one finds, for instance, for the gluonic 
and fermionic contributions to the energy density of QCD with nj massless flavours 



On lattices with small temporal extent these 0{N^'^) corrections lead to large dis- 
tortions of the continuum Stefan-Boltzmann limit 0. In the staggered fermion 
formulation the deviations from the continuum ideal g large as 77% on a 

lattice with temporal extent A^,- = 4 and are still 20 % for = 12. 

As the computational effort for calculating thermodynamic observables, which 
generally have dimension a~^, increases approximately like N^^ it is highly desired 
to be able to extract continuum physics from results of simulations on lattices with 
small temporal extent. In addition the improvement of the fermion action rapidly 
becomes computationally very demanding. In the inversion of the fermion matrix 
the time per iteration is proportional to the number of non-zero entries in a given 
row of the fermion matrix, i.e. the number of neighbours to a given lattice site 
used to discretize the kinetic part of the Lagrangian. In particular in simulations 
with dynamical fermions this gives a very serious constraint and does, at present, 
make it impossible to use, for instance, a complicated fixed point action |Q. In 
this first exploratory study, which aims at an examination of the influence of an 
improved fermion sector on the behaviour of thermodynamic observables at high 
temperature, we therefore have chosen a straightforward improvement scheme for 
the free fermion sector. The minimal extension of the standard staggered fermion 
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discretization scheme which is capable to remove 0(a^) errors in the free fermion 
action involves in addition to the standard one-link term an appropriately weighted, 
straight three-link term. This preserves all the symmetries of the staggered action 
1^ and doubles the computational effort in the fermion sector. 

In the next section we specify the improved action used here and discuss the 
cut-off dependence of thermodynamic observables in the infinite temperature, ideal 
gas limit. Section 3 is devoted to a discussion of the equation of state of four-flavour 
QCD. In section 4 we discuss a calculation of the critical temperature in units of 
the string tension. Finally we present our conclusions in section 5. 



2 Improved actions 



When formulating a discretized version of QCD one has a great deal of freedom in 
choosing a lattice action. Different formulations may differ by subleading powers 
of the lattice cut-off, which vanish in the continuum limit. This has, for instance, 
been used to systematically improve lattice regularized SU{N) gauge theories fol- 
lowing the Symanzik improvement programme 0] • In addition to the elementary 
plaquette term appearing in the standard Wilson formulation of lattice QCD larger 
loops can be added to the action in such a way that the leading O(a^) deviations are 
eliminated from the free fermion action. In our previous studies of the pure gauge 
sector of QCD [|l], Q we have found that already the addition of a planar (1,2) loop 
to the standard Wilson action eliminates the main systematic errors in the high 
temperature limit and even at Tc compares well with a non-perturbatively improved 
(tadpole) action and a fixed point action 0. We therefore use the tree-level 
improved (1,2)-Symanzik action for discretizing the pure gauge sector of QCD, 




In the fermion sector we use an improved action, = ipMip, where a higher order 
difference scheme is used to eliminate the 0{a'^) errors in the discretization of the 
derivatives 9^ appearing in the free fermion action. A three-link term is added 
which preserves all the symmetries of the staggered action . Gauge fields are then 
introduced on the shortest path connecting the staggered fermion and anti-fermion 
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fields. The improved fermion matrix thus reads 

M[U],,= m5,, + r/. ( [t/],, - ^5 [f/],,) (3) 
with Tji denoting the phase factors for staggered fermions and 

The entire action is given by S^[U] = PS^^'"^^ + S[. 

The importance of improved actions for thermodynamic calculations becomes 
evident from an analysis of the high temperature ideal gas limit on lattices of size 
N^Nr- One indeed finds a strong reduction of the cut-off dependence relative to the 
standard Wilson formulation. In order to quantify this we previously had calculated 
the energy density of a free gluon gas ^ which shows explicitly that corrections 
to the continuum Stefan-Boltzmann law only start at 0{N~'^). We give here the 
corresponding results for a free massless fermion gas using the improved fermion 
action for calculations on spatially infinite lattices with temporal extent Nt, 



2^4 



Nr-l 



/2((2no + l)7r/iV.) 



'oCu2(27rp)+4/2((2no + l)vr/iV,) 



' /^(27rpo) 
^°o^2(2vrp)+4/2(27rpo) 



(5) 



where no = 0, 1,..., A^,- — 1 labels the discrete set of Matsubara modes, uj'^{2ti'p) = 
4=J2^=i PiP/j), and the function f{p) gives the momentum dependent terms in the 
free fermion propagator, 



9 1 

fip) = si^b) - ^ sin(3p) . 



(6) 



Results for the free gluon and fermion energy densities are given in Table |1|. We 
note that in the gluon sector the deviation from the continuum limit is below 3% 
already on iV,- = 4 lattices while the fermion action still shows about 27% deviations 
on this size lattice. Still the cut-off dependence is drastically reduced compared to 
the standard staggered fermion formulation, which is given for comparison in the 
last column of Table |l[ 
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Table 1: The energy density of a free gluon (e*-^) and a free fermion (e"^) gas on spa- 
tially infinite lattices with temporal extent A^,- relative to the corresponding values 
in the continuum, e'^/T^ = Svr^/lS and e^/T^ = 7?7,j7r^/20. The last two columns 
give results for the energy density of four-flavour QCD in the infinite temperature 
limit relative to the corresponding continuum result, esslT^ = (e*^ + e^)/T*. 





0(0^) improved 


standard 




e^(iV.)/e^ 


e^(iV.)/e^ 


e{Nr)/esB 


4 


0.986568 


1.269896 


1.191737 


1.694917 


6 


0.997528 


1.005633 


1.003397 


1.659614 


8 


1.000309 


0.963731 


0.973822 


1.424757 


10 


1.000253 


0.986795 


0.990507 


1.252452 


12 


1.000150 


0.996509 


0.997513 


1.156253 



3 QCD equation of state 

We have analyzed the thermodynamics of four-flavour QCD using the improved 
action defined in Eq. ^ We have performed calculations on 16^ x 4 and 16^ lattices 
with quark masses ma = 0.05 and 0.1 using the hybrid Monte Carlo algorithm. The 
inclusion of the contributions from the three-link terms to the equations of motion is 
straightforward. At each value of the gauge coupling, /?, typically 1500 trajectories 
of length At = 0.6 have been analyzed on the finite temperature lattice and 500 
trajectories on the zero temperature lattice. 

The basic observables entering the calculation of thermodynamic quantities are 
the expectation values for the gauge action, (S''-^'^^), and the chiral condensate, 
(xx)) ^-C- the derivatives of the partition function Z = J df/d^d^exp (— 5^) with 
respect to the gauge coupling and the quark mass, respectively. In particular we will 
need the differences between expectation values calculated on the finite temperature 
lattice (16^ x 4) and the zero temperature lattice (16^), 

ixx) = {xx)o - {xx)t , (7) 

where the normalization of the expectation values in both cases has been defined 
per lattice site, z.e. {S^^'^^) = -N-^N-^d\n Z/df3 and (xx) = N-^ N-^ d In Z/dma. 
Our results for these differences are summarized in Figure |I| for both values of the 
quark mass. The curves shown in the Figure are spline interpolations, which have 
been used for the subsequent calculation of the pressure {p) and energy density (e) . 
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Figure 1: Difference of expectation values 
densates calculated on 16^ x 4 and 16^ lat 
ma = 0.05 and 0.1. 




0.5 1 1.5 2 2.5 3 3.5 

of the gluonic action and chiral con- 
ices at two values of the quark mass, 



Like in the pure gauge theory the pressure can be obtained from an integration 



of the difference of gluonic action densities [|1^ 



P_ 

2^4 



/3o 







Nt I d/3'(^(i'2)) 
ho 



(8) 



The calculation of energy density involves in addition also the differences of the 
chiral condensates at zero and non-zero temperature. Furthermore one needs to 
know the cut-off dependence of the two bare couplings, (3 and ma. 



A 



3p 



2^4 



-Nt 



(9) 



with 



Rf-. 



d/3 
din a 



Rr, 



dma 
din a 



(10) 
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In the chiral limit the derivative Rm vanishes and (e — 3p) is again proportional only 
to the /3-function, i?^, as it is the case in the pure gauge sector. The derivatives 
of the bare couplings (3 and ma with respect to the lattice cut-off can be obtained 
from the cut-off dependence of two physical observables. In general we may follow 
Ref. |jl2| and use two meson masses for this, m^r and rup. If we parametrize the 



quark mass dependence of these masses as 

a^ml = h{P)ma , 
amp = f{f3) + g{/3)ma , (11) 

we find for the derivatives 

f 1 - maf f 

= 77TT f,' k',g =ij + Oima) , 

f'l+ma{^-f)f f 

Rm = 2ma(^l-^Rf,^ . (12) 

Here the prime denotes derivatives with respect to (3. In the asymptotic scaling 
regime the ratios ///' as well as h/h' will be given by the leading term in the QCD 
renormalization group equation, i.e. f / f = h/h' = — 25/47r^. The second derivative 
thus becomes asymptotically Rm = ma + 0{{ma)'^) . We find strong deviations from 
these asymptotic relations in the parameter range studied here. This is immediately 
evident from the quark mass dependence of the pion mass. To a first approximation 
the pion slope turns out to be independent of (3. This shows that Rm — 2ma. In 
fact, for both quark mass values and all /5- values studied by us we find this relation 
to be satisfied within 15%. 

Rather than using the rho-meson mass to determine Rp we have used an approach 
which turned out to describe quite well the deviations from asymptotic scaling in the 
pure gauge theory. We define an effective coupling, P^s, in terms of the expectation 
value of the gluonic part of the action, 

12 

^''''^WW) ' ^^^^ 

and use the dependence of on P to calculate the derivative Rj3 with the help 
of the asymptotic two-loop renormalization group equation. This also fixes the 
temperature scale, 

T / B \ -77/625 

^= (^) exp(47rV25(/3efr -/?.)) , (14) 
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We use this approximation for Rf^ also in the calculation of Rm and determine in 
addition the ratio h'/h from the pion slope. 

The critical couplings for both quark masses have been determined from the 
Polyakov-loop susceptibility and the chiral susceptibility. The location of the peaks 
in both susceptibilities coincide within errors. From this we find 



3.57 (3) 
3.49 (3) 



ma 
ma 



0.1 
0.05 



(15) 



which fixes the temperature scale through Eq. IT| 



In Figure ^ we show the results of a calculation of the pressure for both values 
of the quark mass. As can be seen the result is quite insensitive to the value of 
the quark mass. Moreover, we note that the overall structure of the temperature 



dependence of the pressure is quite similar to that of the pure gauge sector |TT 
This also can be seen in Figure 
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where we show the pressure of the SU{3) gauge 



P/T^ 






Psb{N,=4)/T^ - 




'^'''''^^ m/T=0.4 




\\ m/T=0.2 - 




\ quenched 
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Figure 2: Comparison between the pressure of four-flavour QCD on a 16^ x 4 lattice 
for two values of the quark mass and the pressure of the pure SU (3) gauge theory. 
The latter has been rescaled by the appropriate number of degrees of freedom of 
four- flavour QCD (29/8) and a factor 1.19 which takes care of the remaining finite 
cut-off distortion of the ideal gas limit in our four-flavour simulation (see Table 1). 

theory rescaled by an appropriate ratio of the number of degrees of freedom so that 
the high temperature limit coincides with that of four- flavour QCD. 

In Figure ^ we show the difference between the energy density and three times 
the pressure calculated according to Eq. ||. As discussed above, the second term in 
this equation, i.e. the term being proportional to the chiral condensates, will not 
contribute in the chiral limit. We therefore show in Figure |^ also the difference 
((e — 3p)/T^)o defined only in terms of the gluonic part, 

Ao = -N^RpW^ (16) 
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Figure 3: The difference (e — 3p)/T'^ calculated according to Eq. 9 (A) as well as 
Eq. 16 (Ao).Tlie lines are spline interpolations to the data, with the full line for 
m/T = 0.2, the dotted one for m/T = 0.4 and dash-dotted for the quenched results. 

We note that in the high temperature regime, T/Tc>2, the quark mass depen- 
dence, which is visible in A is completely eliminated in Aq. This suggests that at 
least in this temperature regime Aq can be viewed as a reasonable extrapolation of 
(e — 3p)/T'^ to the chiral limit, ma — > 0. Close to Tc this is not that evident. In 
fact, neither in Figure ^ nor in |]b we see any significant quark mass dependence 
for T<1.2Tc. In the case of (e — 3p)/T^ this is due to the cancellation of two effects 
in the contribution from the chiral condensate term. On the one hand the derivative 
factor Rm decreases with decreasing quark mass. On the other hand the difference 
of chiral condensates, {xx)y still increases because the critical coupling gets shifted 
towards smaller values (Eq. p!5| ). As the first effect will ultimately dominate and will 
force the entire contribution of the condensate term to vanish in the chiral limit we 
expect that Aq is already a good approximation to the final result in the chiral limit. 
Of course, this should eventually be checked by performing simulations at smaller 
quark masses. 

Also shown in Figure ^ is the result for the pure gauge theory which again has 
been rescaled with the appropriate number of degrees of freedom (29/8) as it has 
been done for the pressure alone. This too gives support to our extrapolation to 
the chiral limit and suggests that the equation of state in four-flavour QCD indeed 
shows a temperature dependence which is very similar to that of the pure gauge 
theory. 

We use the results for the pressure and the difference e — 3p to extract the energy 
density. This is shown in Figure ^. The energy density does come close to the ideal 
gas limit immediately above Tc. We do observe an overshooting of the ideal gas limit 
close to Tc for the non-zero quark masses considered here. This is a feature not seen 
in the pure gauge sector. The overshooting does, however, seem to disappear when 
we determine the energy density from Aq. As discussed above this may be viewed 



as an extrapolation to the chiral limit (see also [|16 
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Figure 4: Energy density of four-flavour QCD on a 16^ x 4 lattice. The lower set 
of curves shows an extrapolation to the chiral limit which has been obtained by 
ignoring the second term in Eq. 9, i.e. from Eq. 16 (see text). 

4 The critical temperature 



In the pure gauge sector it has been found that the critical temperature for the 
deconfinement transition calculated with the standard Wilson action in units of the 
square root of the string tension, ^/o^, deviates by about 10% from the continuum 
extrapolation at aT^ = 1/4. This cut-off dependence gets strongly reduced already 
with tree level improved actions [13, n]. To our knowledge the ratio Tc/y/a has 



previously been determined for the chiral transition in four-flavour QCD only once 
using the standard staggered formulation with a quark mass m/T = 0.08 at aTc = 
1/8 [|15[. This gave TJ^/^ = 0.39 (3). 

We have analyzed the heavy quark potential on the 16^ lattices at the critical 
couplings given in Eq. |T5|. For this purpose we have generated 300 gauge field 
configurations separated by 10 trajectories for the quark mass ma = 0.05 and 200 
configurations for the larger mass, ma = 0.01. The potential has then been extracted 
using smeared Wilson loops and following the approach used also in the pure gauge 
theory [|1^]. The potentials have been fitted to an ansatz including a linear and a 
Coulomb term, Vqq{R) = Vo + a/R + aR. The string tension in units of the cut-off, 
a = aa^ turns out to be significantly larger than in the SU(3) gauge theory at a 
comparable value of the cut-off, i.e. aTf. = 1/4. The potentials have been fitted in 
the range 3 < i? < 6. From the fits we obtain 



'(7 



0.407 (10) 
0.430 { 8) 



ma 
ma - 



- 0.05 
0.1 



(17) 
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5 Conclusions 



We have analyzed the thermodynamics of four-flavour QCD using an improved 
fermion action. We find that the O(a^) improvement of the free fermion action 
also leads to a strong reduction of the cut-off dependence in the high temperature 
phase. In the ideal gas limit deviations from the continuum Stefan-Boltzmann law 
are still about 20%. A further improvement will thus be necessary in order to reduce 
the cut-off dependence to only a few percent as it is the case in the pure gauge sec- 
tor. In addition one would also like to achieve a reduction of the flavour symmetry 
breaking in the staggered action. This does not seem to be the case for the improved 



action we have used here \l7j. 



For the dependence of the pressure on T/T^ we obtain results which closely 
follow that of the SU{3) gauge theory when rescaled with the number of degrees of 
freedom. The same holds for the energy density after an extrapolation to the chiral 
limit. The critical temperature itself drops to = 0.407(10)^/o" at a quark mass 
m/T = 0.2. 



References 

[1] B. Beinlich, F. Karsch and E. Laermann, Nucl. Phys. B462 (1996) 415. 
[2] A. Papa, Nucl. Phys. B478 (1996) 335. 

[3] B. Beinlich, F. Karsch and A. Peikert, SU(3) Latent Heat and Surface Tension 
from Tree Level and Tadpole Improved Actions, BI-TP 96/24 (|hep-lat/960814lD , 
to appear in Phys. Lett. B. 



[4] W. Bietenholz, R. Brower, S. Chandrasekharan and U. J. Wiese, Progress on 
Perfect Lattice Actions for QCD, |hep-lat /9608068 , to appear in the Proceedings 



of the "International Symposium Lattice '96". 

[5] F. Karsch, J. Engels and H. Satz, Nucl. Phys. B205 [FS5] (1982) 239. 

[6] K. Symanzik, Nucl. Phys. B226 (1983) 187 and Nucl. Phys. B226 (1983) 205. 

[7] P. Weisz, Nucl. Phys. B212 (1983) 1; P. Weisz and R. Wohlert, Nucl. Phys. B236 
(1984) 397 and Nucl. Phys. B247 (1984) 544. 

[8] S. Naik, Nucl. Phys. B316 (1989) 238. 

[9] G.P. Lepage and P.B Mackenzie, Phys. Rev. D48 (1993) 2250. 



11 



[10] J. Engels, J. Fingberg, F. Karsch, D. Miller and M. Weber, Phys. Lett. B 252 
(1990) 625. 

[11] G. Boyd et al, Phys. Rev. Lett. 75 (1995) 4169 and Nucl. Phys. B469 (1996) 
419. 

[12] T. Blum, L. Karkkainen, D. Toussaint and S. Gottheb, Phys. Rev. D51 (1995) 
5153. 

[13] F. Karsch with B. Beinhch, J. Engels, R. Joswig, E. Laermann, A. Peikert 
and B. Petersson, BI-TP 96/34 ( [hep-lat / 9608047|) , QCD Thermodynamics with 
Improved Actions, to appear in the Proceedings of the "International Symposium 
Lattice '96" . 

[14] B. Beinlich, F. Karsch, E. Laermann and A. Peikert, String tension and Ther- 
modynamics from Tree Level and Tadpole Improved Actions, in preparation. 

[15] K.D. Born et al., Phys. Lett. B329 (1994) 325. 

[16] C. Bernard et al.. Thermodynamics for two flavor QCD, [hep-lat /9608026| , to 
appear in the Proceedings of the "International Symposium Lattice '96". 

[17] C. Bernard et al.. Recent MILC spectrum results, [hep-lat /9608102| , to appear 
in the Proceedings of the "International Symposium Lattice '96". 



12 



